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The recent observational evidences for the present accelerated stage of the Universe have stim- 
ulated renewed interest for alternative cosmologies. In general, such models contain an unknown 
negative-pressure dark component that explains the supernova results and reconciles the inflationary 
flatness prediction (fix = 1) and the cosmic microwave background measurements with the dynam- 
ical estimates of the quantity of matter in the Universe (fi m ~ 0.3 ± 0.1). In this paper we study 
some observational consequences of a dark energy candidate, the so-called generalized Chaplygin gas 
which is characterized by an equation of state pc = —A/p c , where A and a are positive constants. 
We investigate the prospects for constraining the equation of state of this dark energy component by 
combining Chandra observations of the X-ray luminosity of galaxy clusters, independent measure- 
ments of the baryonic matter density, the latest measurements of the Hubble parameter as given by 
the HST Key Project and data of the Supernova Cosmology Project. We show that very stringent 
constraints on the model parameters can be obtained from this combination of observational data. 

PACS numbers: 98.80.Es; 95.35.+d; 98.62.Sb 



I. INTRODUCTION 



One of the most important goals of current cosmolog- 
ical studies is to unveil the nature of the so-called dark 
energy or quintessence, the exotic negative-pressure com- 
ponent responsible for the accelerating expansion of our 
Universe. Over the last years, a number of candidates 
for this dark energy have been proposed in the litera- 
ture 0, with the vacuum energy density (or cosmolog- 
ical constant) and a dynamical scalar field 0, Q appar- 
ently constituting the most plausible explanations. From 
the observational viewpoint, these two class of models 
are currently considered our best description of the ob- 
served Universe whereas from the theoretical viewpoint 
they usually face fine-tunning problems, notably the so- 
called cosmological constant problem 4\ as well as the 
cosmic coincidence problem, i.e., the question of explain- 
ing why the vacuum energy or the scalar field dominates 
the Universe only very recently. The later problem hap- 
pens even for tracker versions of scalar field models in 
which the evolution of the dark energy density is fairly 
independent of initial conditions 0, U 

Among the many dark energy candidates, a recent and 
very interesting proposal has been suggested by Kamen- 
shchik et al. _m and developed by Bilic et al. Q and 
Bento et al. Q- It refers to the so-called Chaplygin gas 
(C), an exotic fluid whose equation of state is given by 
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with a = 1 and A a positive constant. Actually, the 
above equation for a^l generalizes the original Chap- 
lygin equation of state proposed in Ref. Q whereas for 
a = 0, the model behaves like scenarios with cold dark 
matter plus a cosmological constant (ACDM). 

In the context of the Friedman-Robertson- Walker 
(FRW) cosmologies, whether one inserts Eq. (1) into 
the energy conservation law {u^Tf^ = 0), the following 
expression for the energy density is immediately obtained 



Pc 



or, equivalently, 
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where pc a is the current energy density (from now on 
a subscript o means present day quantities, and C de- 
notes either the Chaplygin gas or its generalized ver- 
sion). The function R(t) is the cosmic scale factor, 
B = /3p +Q — A is a constant and A s = Aj pQ a is a quan- 
tity related to the present day Chaplygin adiabatic sound 
speed (Ug = aA/p c i+ a ), As can be seen from the above 
equations, the C-gas interpolates between non-relativistic 
matter (pc(R —> 0) ~ VB/R 3 ) and negative-pressure 
dark component regimes (pc(R — » oo) ~ y/A). This 
particular behavior of the Chaplygin gas inspired some 
authors to propose a unified scheme for the cosmologi- 
cal "dark sector" 0, H, @ , an interesting idea which has 
also been considered in many different contexts [ljj (see, 
however, pdjb 
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In the theoretical front, a connection between the 
Chaplygin equation of state and string theory had long 
been identified by Bordemann and Hoppe | 12| and Hoppe 
[l3| (see also Q for a detailed review). As explained in 
such references, a Chaplygin gas-type equation of state is 
associated with a parametric description of the invariant 
Nambu-Goto d-brane action in a d + 2 spacctimc. In the 
light-cone parameterization, such an action reduces itself 
to the action of a Newtonian fluid which obeys Eq. (1) 
with a = 1, with the C-gas corresponding effectively to 
a cZ-branes gas in a (d + 2)-dimensional spacetime. 

Another interesting connection is related to recent at- 
tempts of describing the dark energy component through 
the original Chaplygin gas or its generalized version. 
Such a possibility has provoked a growing interest for 
exploring the observational consequences of this fluid in 
the cosmological context. For example, Fabris et al. 0] 
analyzed some consequences of such scenario using type 
la supernovae data (SNe la). Their results indicate that 
a Universe completely dominated by the Chaplygin gas 
is favored when compared with ACDM models. Re- 
cently, Avelino et al. ^(| used a larger sample of SNe 
la and the shape of the matter power spectrum to show 
that such data restrict the model to a behaviour that 
closely matches that of a ACDM models while Bento et 
al. 0, 0| showed that the location of the CMB peaks 
imposes tight constraints on the free parameters of the 
model. More recently, Dev, Alcaniz & Jain 0] and Al- 
caniz, Jain & Dev [!j investigated the constraints on the 
C-gas equation of state from strong lensing statistics and 
high-z age estimates, respectively, while Silva & Berto- 
lami [2(j studied the use of future SNAP data together 
with the result of searches for strong gravitational lenses 
in future large quasar surveys to constrain C-gas mod- 
els. Makler et al. |2l| also showed that such models are 
consistent with current SNe la data for a broad range 
of parameters. The trajectories of statefinder parame- 
ters in this class of scenarios were studied in Ref. 
[23| while constraints involving Cosmic Microwave Back- 
ground (CMB) data have also been extensively discussed 
by many authors d E El • 

In this work we study the possibility of constraining 
the generalized Chaplygin equation of state from X-ray 
luminosity of galaxy clusters. With basis on measure- 
ments of the mean baryonic mass fraction in clusters as 
a function of redshift, we consider the method originally 
proposed by Sasaki l2t| and Pen [2?|], and further mod- 
ified by Allen et al. |28l |29| who analyzed the X-ray 
observations in some relaxed lensing clusters observed 
with Chandra in the redshift interval 0.1 < z < 0.5 (see 
also H3). By inferring the corresponding gas mass frac- 
tion, Allen and collaborators placed observational limits 
on the total matter density parameter, fl m , as well as 
on the density parameter associated to the vacuum 
energy density. More recently, a similar analysis has also 
been applied to conventional quintessence models with 
an equation of state p x = ujp x by Lima et al. |31|. 

The paper is organized as follows. In Sec. II we present 



the field equations and distance formulas necessary to 
our analysis. In Sec. Ill the corresponding limits on C- 
gas models from X-ray luminosity of galaxy clusters are 
derived. We also examine the limits from a statistical 
combination between X-ray data and recent SNe la ob- 
servations. Finally, in Sec. IV, we finish the paper by 
summarizing the main results and comparing our con- 
straints with others derived from independent analyses. 



II. THE CHAPLYGIN GAS MODEL 

The FRW equation for a spatially flat, homogeneous, 
and isotropic scenarios driven by nonrelativistic matter 
and a separately conserved C-gas component reads 



(|) 2 = Hl{n m {^f 
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where an overdot denotes time derivative, H = 
100/iKm.s~ 1 .Mpc~ 1 is the present value of the Hubble 
parameter, f2 m is the matter density parameter, and the 
dependence of the C-gas energy density with the scale 
factor comes from Eq. (3). 

The comoving distance r\(z) to a light source located 
at r = 7'i and t — t\ and observed at r = and t = t Q is 
given by 



r\{z) = 
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dx 
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where x' = R(t)/R = (1 + z) 1 is a convenient 
integration variable and the dimensionless function 
T{x, Q m , A s , a) is given by 
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Now, in order to derive the constraints from X-ray gas 
mass fraction on the C-gas, let us consider the concept 
of angular diameter distance, Da{z). Such a quantity is 
defined as the ratio of the source diameter to its angular 
diameter, i.e., 



D A = - = R(h) ri - (1 + z)- l R n(z), 
which provides, when combined with Eq. (5), 
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As one may check, for A s = and a — 1 the above ex- 
pressions reduce to the standard cold dark matter model 
(SCDM). In this case, the angular diameter distance can 
be written as 
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FIG. 1: The model function /^° d as a function of the redshift 
for selected values of A a and a and fixed values of Q m = 0.3, 
Q, b h 2 = 0.0205 and h = 0.72. 



III. LIMITS FROM X-RAY GAS MASS 
FRACTION 

Following Allen et al. jH El and Lima et al. [3l| . 
we consider the Chandra data consisting of six clusters 
distributed over the redshift interval 0.1 < z < 0.5. The 
data are constituted of regular, relatively relaxed sys- 
tems for which independent confirmation of the matter 
density parameter results is available from gravitational 
lensing studies. The X-ray gas mass fraction (/ gas ) val- 
ues were determined for a canonical radius ^2500, which 
is defined as the radius within which the mean mass den- 
sity is 2500 times the critical density of the Universe at 
the redshift of the cluster. In order to generate the data 
set the SCDM model with H a = SOKm.s^.Mpc" 1 was 
used as the default cosmology (see (3 for details). 

By assuming that the baryonic mass fraction in galaxy 
clusters provides a fair sample of the distribution of 
baryons at large scale (see, for instance, (H) and that 
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|2q. the model function is defined as [28 
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where the bias factor b ~ 0.93 |3j( is a parameter moti- 
vated by gas dynamical simulations that takes into ac- 
count the fact that the baryon fraction in clusters is 
slightly depressed with respect to the Universe as a whole 
[23. The term (2ft,) 3 / 2 represents the change in the 
Hubble parameter between the default cosmology and 
quintessence scenarios and the ratio D^ BM (z;) / D^(zi) 
accounts for deviations in the geometry of the universe 
from the default cosmology (SCDM model). In Fig. 1 we 
show the behavior of /™° d as a function of the redshift 



FIG. 2: a) The Ax 2 contours for the a — A B plane according to 
the X-ray data discussed in the text. The contours correspond 
to 68%, 95% and 99% confidence levels. The value of the 
matter density parameter has been fixed at fi m = 0.3. b) 
fi m — A s plane for the original C-gas model (a = 1). At 
95.4% we find f2 m = 0.3 ± 0.02 while the entire range of A 3 
is allowed. Note that the X-ray data constrain tightly the 
matter density parameter. 



for some selected values of A s and a having the values of 
f2b and h fixed. For the sake of comparison, the current 
favored cosmological model, namely, a flat scenario with 
70% of the critical energy density dominated by a cos- 
mological constant (ACDM) is also shown. In order to 
have bidimensional plots we fix the value of f2 m as sug- 
gested by dynamical estimates [35|, i.e., 0.3 in Fig. 1 as 
well as in all statistical analyses involving the generalized 
Chaplygin gas. However, in the case of a conventional C- 
gas (a — 1), the density parameter fl m has also been 
considered a free parameter to be adjusted by the data. 

The cosmological parameters A s and a are determined 
by using a % 2 minimization with the priors fl^h 2 — 
0.0205 ± 0.0018 and h = 0.72 ± 0.08 [jtfl for the 
range of A s and a spanning the interval [0,1] in steps of 
0.02, 
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where 07 aa ( are the symmetric root-mean-square errors 
for the SCDM data. The 68.3% and 95.4% confidence 
levels are defined by the conventional two-parameters \ 2 
levels 2.30 and 6.17, respectively. 

In Fig. 2a we show contours of constant likelihood 
(68%, 95% and 99%) in the parameter space a — A s for 
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FIG. 3: a) The likelihood contours in the a — A s plane for 
the joint X-ray + SNe la analysis described in the text. The 
contours correspond to 68%, 95% and 99% confidence levels, 
b) The Q m — A s plane for the joint X-ray + SNe la analysis. 
The best fit values are located at A s = 0.98 and fi m = 0.3 At 
95.4% we find A s > 0.84 and 0.273 < fi m < 0.329. 



the X-ray data discussed earlier. From the above equa- 
tion we find that the best fit model occurs for A s = 1 
which, according to Eq. (4), is independent of the index 
a and equivalent to a ACDM universe. Such model cor- 
responds to a accelerating scenario with the deceleration 
parameter q a = —0.55 |46j . From this figure, we also see 
that both A s and a are quite insensitive to these data 
and that, at 95.4% c.L, one can limit the parameter A s 
to be > 0.52. Figure 2b shows the plane J7 m — A s for the 
conventional C-gas (a = 1). As one should expect from 
different analyses [23, EI > the matter density parameter 
is very well constrained by this data set while the param- 
eter A s keeps quite insensitive to it. The best fit occurs 
for models lying in the interval A s = [0, 1] and f2 ro = 0.3. 
At 95.4% c.L, we find 0.268 < ft m < 0.379. For a X-ray 
analysis where the Cg plays the role of a unified model 
for dark matter /energy, see |3Sj . 

A. Joint analysis with SNe la 

By combining the X-ray and SNe la data sets, more 
stringent constraints on the cosmological parameters f2 m 
and A s are obtained. As it was shown elsewhere, the 
parameter a is highly insensitive to the SNe la data. 
To perform such analysis, we follow the conventional 
magnitude-redshift test (see, for example, jj^) and use 
the SNe la data set that corresponds to the primary fit 
C of Perlmutter et al. ^3 together with the highest su- 
pernova observed so far, i.e, the 1997ff at z = 1.755 and 
effective magnitude m eff = 26.02 ± 0.34 [4l| and two 



FIG. 4: The extended A a — a plane for the joint X-ray + 
SNe la analysis. Although completely modified in comparison 
with the previous ones, the best fit values for this extended 
analysis (A s — 1.02 and a = 0.45) are still in agreement with 
the causality imposed by the sound/light velocity (A s < 1/a). 
The dashed hyperbola corresponds to the limit condition vl = 
1. 



newly discovered SNe la, namely, SN 2002dc at z = 0.475 
and m eff = 22.73 ± 0.23 and SN 2002dd at z = 0.95 and 
m e ff = 24.68 ±0.2 [Z^. Figures 3a, 3b and 4 show 
the results of our analysis. In Fig. 3a we display con- 
tours of the combined likelihood analysis for the para- 
metric space A s — a. In comparison with Fig. 2a we see 
that the available parameter space is reasonably modi- 
fied with the value of A s constrained to be > 0.73 at 
95.4% c.L and the entire interval of a = [0, 1] allowed. 
The best fit model occurs for values of A s — 0.98 and 
a = 0.93 with Xmin — 61-38 and v = 61 degrees of free- 
dom (Xmin/ l/ — 1-0)- The most restrictive limits from 
this joint analysis are obtained for the original version of 
C-gas (a = 1). In this case, the plane f2 m — A s (Fig. 
3b) is tightly constrained with the best fit values lo- 
cated at A s = 0.98 and tt m = 0.3 with xLiJ y - L0 - 
At 95.4% this analysis also provides A s > 0.84 and 
0.273 < Q m < 0.329. Note that the contours a - A s 
(Fig. 3a) and il m — A s (Fig. 2b) are almost orthogonal, 
thereby explaning the shape of the fl m — A s plane ap- 
pearing in Fig. 3b. We also observe that by extending 
the a — A s plane to the interval [0,2] the new best fit 
values (A s = 1.02 and a — 0.45), although completely 
modified in comparison with the previous ones, are still 
in agreement with the causality (A s < 1/cc) imposed by 
the fact that the adiabatic sound speed vl = dp /dp in the 
medium must be lesser than or equal to the light velocity 
(see Eq. 1). Some basic results of the above analysis are 
displayed in Fig. 4. 
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IV. DISCUSSION AND CONCLUSIONS 

Alternative cosmologies with a quintessence compo- 
nent (dark energy) may provide an explanation for the 
present accelerated stage of the universe as suggested by 
the SNe la results. In this work we have focused our at- 
tention on a possible dark energy candidate, the so-called 
Chaplygin gas. The equation of state of this dark energy 
component has been constrained by combining Chandra 
observations of the X-ray luminosity of galaxy clusters 
and independent measurements of the Hubble parame- 
ter and of the baryonic matter density as well as from 
a statistical combination between X-ray data and recent 
SNe la observations. We have shown that stringent con- 
straints on the free parameters of the model, namely A s , 
a and f2 m , can be obtained from this combination of ob- 
servational data. 

It is also interesting to compare the results derived 
here with another independent analyses. For example, 
using only SNe la data, Fabris et al. ^5| found A s = 
0.93^Q2o for the original C-gas model (a = 1) with the 
matter density parameter constrained by the interval < 
< 0.35. The same analysis for f2 m = fib = 0.04 
(in which the C-gas plays the role of both dark matter 
and dark energy) provides A s = 0.871q ]J. These values 
agree at some level with the ones obtained from statistics 
of gravitational lensing (SGL), i.e., A s > 0.72 and 
age estimates of high-z galaxies (OHRG's), A s > 0.85 
- A s > 0.99 for the interval Ct m = 0.2 — 0.4 with lower 
values of A s corresponding to lower fl m [t| . The original 



Chaplygin gas model, however, seems to be incompatible 
with the localization of the acoustic peak of CMB as 
given by WMAP [H and BOOMERANG Q data. For 
the case of a generalized component, the same analysis 
shows that for intermediary values of the spectral til t n, ? , 
the C-gas model is favored by this data set if a ~ 0.2 fisl 
A similar analysis for BOOMERANG and Arqueops |45l| 
data implies 0.57 < A s < 0.91 for a < 1 @ whereas 
an investigation involving WMAP and SNe la data sets 
restricts a to be < a < 0.2 p5|. 

It should be stressed that our results are in line with 
the above quoted independent studies. In particular, 
even considering that the parameter A s is quite insen- 
sitive to the X-ray data alone, the matter density pa- 
rameter is very well constrained. This result is also in 
agreement with the limits derived by Allen et al. 29] for 
ACDM models. In addition, as shown in Figure 3b, by 
combining the X-ray and SNe la data sets, more stringent 
constraints on the parameters Q m and A s are readily ob- 
tained. From the above analyses we also note that the a 
parameter is more strongly restricted if causality require- 
ments (Vg < 1) are imposed (see Figure 4). However, it 
seems that a even better method to place limits on such 
a parameter is through the physics of the perturbations, 
i.e., CMB and LSS data (see, e.g., (HIH). 
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